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SUBSTITUTE SPECIFICATION (Clean Copy) 

METHODS OF ESTIMATING PRECIPITATION CHARACTERISTICS 

Related Application 

[0001] This is a §371 of International Application No. PCT/FR2004/002692, with an inter- 
national filing date of October 21, 2004 (WO 2005/040853 Al, pubhshed May 6, 2005), which 
is based on French Patent Application No. 03/12299, filed October 21, 2003. 

Technical Field 

[0002] This invention relates to methods of estimating precipitation characteristics and, in 
one particular aspect, precipitation rate for solid precipitation. 

Background 

[0003] It is known that the characteristics of a radar image can be used to estimate rainfall 
characteristics. In particular, EP 1 049 944 describes a technique for estimating rainfall using a 
radar. In that patent, provision is made for the following steps: 

using a dual-polarized radar to measure the differential phase shift (Odp) and the 

apparent reflectivity Z in at least one of the H (horizontal) or V (vertical) polarizations 

and over a given range (ro, ri) of path radii relative to the radar; 

determining an estimate of the value No* representative of the size distribution of 

the rain drops, on the basis of the differential phase shift difference over the range ro to 

ri, and on the basis of an integral of a fimction of the apparent reflectivity Z, over the 

range (ro, ri); and 



deducing the value of the precipitation rate at a point on the basis of No* and of 
the apparent reflectivity at the point. 
[0004] A method of estimating rainfall is also known from WO 03/007016, which describes 
a method of estimating a precipitation rate by means of a dual-polarized radar, that method being 
characterized by the following steps: 

using a dual-polarized radar to measure the differential phase shift Odp and atten- 
uated reflectivity Z in at least one of the H or V polarizations and over a given range (rl, 
rO) of path radii r relative to the radar; 

determining an estimate of the value K(rO) of the specific attenuation at rO on the 
basis of the profile of the attenuated reflectivity measured in this way, and on the basis of 
the differential phase shift difference over the range rO to rl; determining an estimate 
K(r) of the specific attenuation at r as a function of the attenuation K(rO) determined in 
this way and of the attenuated reflectivity profile Z(r); and determining the precipitation 
rate R(r) once K(r) is known. 
[0005] Those various solutions make it possible to characterize liquid precipitation and esti- 
mate precisely rainfall rate (in millimeters per hour (mm/h), but they do not make it possible to 
estimate the characteristics of solid precipitation such as snow. 

Summary 

[0006] This invention relates to a method of estimating precipitation characteristics including 
acquiring a radar image including at least a vertical plane of a precipitation zone, processing a 
vertical profile to generate digital signals representative of reflectivity in vertical direction z, 
integrating the signals representative of reflectivity by assimilation of a reflectivity vertical pro- 
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file in an aggregation model to generate a signal representative of the profile in the vertical plane 
of a mean particle diameter weighted by mass of each particle, and determining concentration of 
the solid particles on the basis of signals previously determined. 

Brief Description of the Drawings 

[0007] Selected aspects of the invention will be better understood on reading the following 
description, given with reference to a non-limiting example and selected graphic representations 
below: 

Fig. 1 is a graph of a vertical profile of Z to be inverted (in this example, the isotherm 
0°C is at ground level); 

Fig. 2 is a graph of a profile of Dm resulting fi-om Z being inverted by the aggregation 
model compared with a conventional estimator; 

Fig. 3 is a graph of profiles of NO and of nT resulting fi*om Z being inverted by the aggre- 
gation model compared with a conventional hypothesis and observations; 
Fig. 4 is a graph of the vertical profile of Z to be inverted; 

Fig. 5 is a graph of the profile of Dm resulting firom Z being inverted by the aggregation 
model compared with the conventional estimator; 

Fig. 6 is a graph of the profiles of NO and nT resulting firom the Z being inverted by the 
aggregation model compared with a conventional hypothesis and observations; 
Fig. 7 is a graph of the profile of R resuUing fi*om Z being inverted by the aggregation 
model compared with a conventional estimator; 

Fig. 8 is a graph of the sensitivity of the retrieval of the equivalent precipitation rate 
(mm/h) to the particle density law p(D) oc D-y; and 
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Fig. 9 is a graph of the sensitivity of the retrieval of the vertical Doppler velocity to the 
particle density law p(D) oc D-y. 

Detailed Description 

[0008] This disclosure provides methods of estimating precipitation characteristics and, in a 
particular selected aspect, a precipitation rate for solid precipitation, the method comprising an 
acquisition and processing step comprising acquiring a radar image including at least a vertical 
plane of a precipitation zone, and in processing a vertical profile to deliver digital signals repre- 
sentative of reflectivity in the vertical direction h, an integration step comprising integrating sig- 
nals representative of reflectivity to deliver a signal representative of the profile in the vertical 
plane of a mean particle diameter weighted by the mass of each particle, and a determination step 
comprising determining the concentration of the solid particles on the basis of the signals com- 
puted in the preceding steps. 

[0009] Reference will be made to the drawings below as appropriate. Steps in accordance 
with selected aspects of the method may be seen in one sense by reference to Figs. 1-3. 
[0010] In a preferred aspect, the integration step comprises determining the variable Z(h) of 
the radar observable in nraiVm^ as a function of the altitude h on the basis of the radar image, 
and determining the mean diameter Dm(h) of the particles by solving the following equation: 

fDm b-5 Tll/Z ^ 

= -0.25keffaD 10"*^Z+ I Id^, (2) 

fh l6Z/h j 

where: 

Z is the radar observable to be inverted in mm m' ; 
Dm is in meters (m); 
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a and b are coefficients specific to particles of the "aggregate" type; for example, the co- 
efficient a is equal to 35184 and the coefficient b is equal to 3.16; 

^e^is the coefficient of effectiveness of the aggregation process to be adjusted, said coef- 
ficient keff being equal to 0.3. 

[0011] Integration of (2) uses an integration boundary condition. Advantageously, the inte- 
gration boundary condition is determined so that the value Dm(h) at the top of the cloud cor- 
responds to the predetermined value for the total number of particles at the top of the cloud. 
[0012] In an advantageous implementation, the profile of the total number of particles nt(h) 
is determined by the following equation: 

nT(h) = x.Z(h)/Dn,(hf 
where x is equal to 25.4 10"^^. 

[0013] In another implementation, the meteorological parameter No(h) is determined by the 
following equation: 

No(h) = y.Z(h)/D„,(h)' 
where y is equal to 102 10 *^. 

[0014] In another aspect, the meteorological parameter corresponding to the profile of the ice 
water content IWC(hXin g/m^) is determined by the following equation: 

IWC(h) = wZ(h)/D4h)^ 
where w is equal to 1.25 10" . 

[0015] In yet another aspect, the meteorological parameter corresponding to the profile of the 
solid precipitation rate R(h) (mm/h equivalent melted) is determined by the following equation: 

R(h) = r.Z(h)/D„,(h)'-'' 
wherein r is equal to 4.698 10"^^. 
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[0016] The method applies to stratiform precipitation. It considers that ice-forming nuclei 
are activated at highly negative temperatures, i.e., at the top of the cloud. Ice crystals formed at 
high altitude settle and grow as they fall, either by sublimation of the ambient saturating water 
vapor, or by collection and freezing of supermelted cloud water droplets, or by random aggre- 
gation of their collisions with other ice crystals. Of the three growth processes, only aggregation 
changes the ice particle concentration. The characterization method is based on a simplified de- 
scription of the aggregation mechanism. We use a class of "profiler" methods that inverts the 
vertical profile of reflectivity measured in the ice to extract therefrom the vertical profile of the 
solid precipitation rate. 
[0017] Selected steps include: 

1 - Particle size distribution expressed in "equivalent melted diameter" is assumed to be 
exponential, i.e.: 

N(D)=Noexp(-4D/Dni) (1) 
where N(D) is the concentration of particles per cubic meter (m^) and per diameter range, 
and No and Dm are the two parameters that characterize distribution. 

2 - The top hmax and the base hmin of the layer of soHd precipitation are determined; 

a. hmax is the maximum altitude of the measured reflectivity profile Z(h). 

b. hmin is either the altitude of the isotherm C'C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 

3 - The profile of the parameter Dm(h) in the range hmax to hmin is then determined by re- 
solving the following differential equation: 

= -0.25keffaD 10'^^Z+ I Id^ (2) 

/h l6Z/h j 

where: 
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- Z is the radar observable to be inverted in mm m , 

- Dm is in meters (m); 

- a and b are coefficients specific to particles of the "aggregate" type, equal re- 
spectively to 35184 and to 3.16 on the basis of the observations of J.D. LocateUi 
and P.V. Hobbs, Fall speeds and masses of solid precipitation particles^ J. Geo- 
phys. Res., 79, 2185 - 2197 (1974), the subject matter of which is incorporated 
herein by reference; 

- keff is the coefficient of effectiveness of the aggregation process to be adjusted 
(the value keff== 0.3 seems correct). 

4 - The integration of (2) takes place fi'om the top, where the boundary condition is ex- 
pressed by fixing the total number of particles «r (or the number of ice- forming nuclei 
activated at the top of the cloud). It is possible to take «r ^maj = 10^ m""', which makes it 
possible to express the boundary condition Dm(hmax) as: 

Dm(HmJ^25,4 10'''(Z(hmJnj(hmjf' (3) 

5 - Once the profile Dni(h) from hmax to hmin has been determined, the profiles of the other 
parameters of interest are computed by the following expressions: 

a. Profile of A^o: No(h)=102.10-^^Z(h)/Dm(h)^ 

b. Profile of the total number of particles nT(h) (in m"^): NT(h)=102.10"^^Z(h) 
/D^h)' 

c. Profile of the ice water content IWC(h) (in g/m^): 

IWC(h)=l .25. 1 0-*^Z(h)/Dm(h)3 
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d. Profile of the solid precipitation rate R(h) (mm/h equivalent melted). By using 
the terminal fall velocity determined by Locatelli and Hobbs for aggregates: 
(vr=107.6 D^'^^ (D in m)), R(h) is expressed by: 
R(h)=4.698.10'*^Z(h)/Dm(h)^'^^ 

[0018] The following description relates to another selected aspect of a method of processing 

for determining the precipitation rate. 

[0019] For rainfall, the algorithm used is the ZPHI algorithm which is the subject of FR 
98/00714 and FR 01/09206, the subject matter of which is incorporated by reference. 
[0020] For soUd precipitation, the following description applies. 

[0021] The method for estimating the precipitation rate for solid precipitation includes a 
class of "profiler" algorithms that inverts the vertical profile of reflectivity measured in the ice to 
extract therefi-om the vertical profile of the solid precipitation rate. 

[0022] The method appUes to stratiform precipitation. It considers that the ice-forming nuc- 
lei are activated only at highly negative temperatures, i.e., at the top of the cloud. The ice crys- 
tals formed at high altitude settle and grow as they fall, either by sublimation of the ambient sat- 
urating water vapor, or by collection and fi-eezing of supermelted cloud water droplets, or else by 
random aggregation of their collisions with other ice crystals. Of the three growth processes, 
only aggregation changes the ice particle concentration. The algorithm is based essentially on a 
simplified description of the aggregation mechanism. Selected steps are as follows: 

Particle size distribution expressed in "equivalent melted diameter" is assumed to 

be exponential, i.e.: 

N(D)=NoQxp(-4D/DJ (1) 

where D is the equivalent melted diameter of the ice particle; 
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N(D) is the concentration of particles per and per diameter range; and 
No and Dm are the two parameters that characterize distribution. 

The top hmax and the base hmin of the layer of solid precipitation are determined; 
hmax is the maximum altitude of the measured reflectivity profile Z(h). 
hmin is either the aUitude of the isotherm 0°C if the ground-level temperature is 
positive, or it is ground level if the ground-level temperature is negative. 
[0023] The profile of the parameter Dm(h) in the range hmax to hmin is then determined by re- 
solving the following differential equation, whose solution can be determined analytically: 

aOm 4^ kefTJ(b,d)a b-(a-2) 1 5Z 



= -24 ZD + Dm (2) 

ah r(4+d) a(a-l) (a-l)Zah 



where: 



foo fx r 1 

J(b4)= I exp(-4x)x'' I Ix^'-x'*' |exp(-4x')dx'dx 

Jo Jo L J 

describes the collision fi-equency; 

Z is the radar observable to be inverted 
Dm is in m; 

- a, c, and d are coefficients dependent on the density law F(D) for ice particles, 
which law is assumed to vary at i)'^ (where 7 can take values in the range 0.25 to 1.1 
depending on the type of particles). These coefficients are drawn from Mitchell's 
Theory (Journal of Atmospheric Sciences, 53, 12, 1996), the subject matter of which 
is incorporated by reference, for representing the following power laws: 

v,(D) = cD"^ (where V, is the terminal fall velocity of the particle of equivalent 

melted diameter D): 
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A = aD (where A is the effective cross-section of the particle of equiv- 
alent melted diameter D); 
With coefficients depending on fi"equency and as a fimction of (D), representing 
by a power law the relationship between reflectivity (measured parameter), the concentra- 
tion of the particles, and their mean equivalent diameter: 
Z = aNoDj 

keff is the coefficient of effectiveness of the aggregation process to be adjusted 
(the value keff = 0.3 seems correct immediately above the isotherm O^C). 
[0024] The integration of (2) takes place fi-om the top, where the boundary condition is ex- 
pressed by fixing the total number of particles nr (or the number of ice-forming nuclei activated 
at the top of the cloud). It is possible to take ht (hmax) = 10^ m'^, which makes it possible to ex- 
press the boundary condition Dni(hmax) as: 

1 

Dm(hmJ= I I (3) 

[0025] Once the profile Dm(h) fi-om hmax to hmin has been determined, the profiles of the other 
parameters of interest are computed by the following expressions: 
Profile oiNo: 
Z(h) 

No(h)= : — 

Profile of the total number of particles n-rCh) (in m'^): ni(h)=02SNo(h)Dm(h) 
Profile of the ice water content IWC(h) (in g/m^): 

IWC(h)'={\ .2272) 1 0'*No(h)D„,(h/ 
Profile of the solid precipitation rate R(h) (mm/h equivalent melted): 
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r(4+d) 

7?^;=(1.885)10^ c No(h)Dn,(h/'^^^^ 

4(4+d) 

Vertical profile of the Doppler velocity VoCh) (in meters per second (m/s)), given 
by a power law dependent on density _(D): 

VD(h)=pDm(hf 

[0026] Because of the dependency of the coefficients relative to the density law F(D), the re- 
trieval of the precipitation rate depends critically on the parameter 7, as shown by Fig. 7. The 
density law thus constitutes a key parameter. 

[0027] By operating in vertical firing, it is possible to measure the Doppler velocity VoCh), 
and compare it with VoCh) computed by the algorithm. The Doppler velocity profile is discrim- 
inating, as shown by Fig. 8. By successive adjustments, it is thus possible to determine the par- 
ameter 7 of the law F(D) for which consistency between VoCh) as computed and VoCh) as mea- 
sured is achieved. Fig. 4 shows an example of a profile inverted by the algorithm (in that exam- 
ple, the isotherm 0°C is at ground level). Fig. 5 shows curves resulting firom application of the 
aggregation model and a conventional model, making it possible for the profile of Dm resulting 
fi"om the inversion of Z by the aggregation model to be compared with the conventional esti- 
mator. 

[0028] Fig. 6 shows curves of the profiles as a fiinction of altitude (m) above the isotherm 
0°C. It enables the profiles of No and ni resulting fi-om Z being inverted by the aggregation ^ 
model to be compared with the conventional hypothesis and observations. Fig. 7 shows the pro- 
file of R resulting fi-om Z being inverted by the aggregation model compared with the conven- 
tional estimator. Fig. 8 shows the sensitivity of the retrieval of the equivalent participation rate 
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(mm/h) to the particle density law p(D) oc D"^. Fig. 9 shows the sensitivity of the retrieval of the 
vertical Doppler velocity to the particle density law p(D) cx; D'^. 
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